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Abstract. We introduce a new family of probability distributions on the set of pure states 
of a finite dimensional quantum system. Without any a priori assumptions, the most natural 
measure on the set of pure state is the uniform (or Haar) measure. Our family of measures 
is indexed by a time parameter t and interpolates between a deterministic measure (t = 0) 
and the uniform measure (t — oo). The measures are constructed using a Brownian motion 
on the unitary group Un- Remarkably, these measures have a Un—i invariance, whereas the 
usual uniform measure has a Un invariance. We compute several averages with respect to 
these measures using as a tool the Laplace transform of the coordinates. 



1. Introduction 

Defining models of randomness for quantum objects has become a central problem in quan- 
tum information theory which has found many interesting and novel applications. Probabil- 
ity measures on the set of quantum states have been investigated thoroughly in recent years 
[17, 14, 2] both in the physical and the mathematical literature. Random quantum channels 
have also been a subject of interest [10, 5, 6]. In particular, ensembles of quantum channels 
are the central idea behind the recent breakthroughs in the additivity conjecture [9]. Random 
quantum states can arise in two different ways. First, they describe states of open systems 
which are subjected to random interaction with an (unknown) environment. This aspect has 
been the starting point of the so-called induced measures which describe finite dimensional 
systems in interaction with a usually larger, but finite dimensional environment. Statistical 
ensembles of quantum states can also be used to study physical properties of generic states, 
such as entanglement, purity or other physically relevant quantities. 

Defining a model of randomness for quantum states amounts to specifying a probability 
measure on the set of density matrices on the corresponding Hilbert space. When one considers 
only pure states (rank one density matrices), it turns out that there exists a unique natural 
candidate for such a probability measure: the Lebesgue measure on the unit sphere of the 
underlying Hilbert space, ir^. This measure on the set unit vectors, called the Fubini-Study 
or the uniform distribution, is canonical in the sense that it is invariant under changes of 
bases: an element ijj of this ensemble has the same distribution as one of its rotations Ui/j, 
for any given unitary matrix U £ £/jv(C). This invariance property, which characterizes the 
uniform distribution, justifies its use when no information about the internal structure of the 
system is known. 

Very recently, new ensembles of pure states have been constructed, in order to take into 
account any available a priori information on the system. In [7], the authors introduce an 
ensemble of states for quantum multipartite systems. Given a graph which encodes the en- 
tanglement between the different parties, a probability measure on the pure states of the total 
system is constructed. The ensemble of graph states is different from the uniform ensemble 
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because it contains the structure information about the initial entanglement present in the 
system. 

In the present work, we are going to generalize the uniform ensemble in another direction, 
by removing the rotation invariance condition. Pure state ensembles on a single-partite system 
with no full invariance property exhibit preferred states which have a larger probability than 
other states. The model introduced in the present work has a symmetry group of smaller 
dimension Un-i(C) and this feature makes is suitable for modeling systems on which some 
partial information is available. Given a fixed state ip of a quantum system, we shall introduce 
a one parameter family of probability measures 7r t indexed by a real parameter t ^ 0. The 
parameter t can be interpreted as a time parameter is such a way that for t = 0, the measure 
is deterministic, being supported on the state ip, and in the limit t — >■ oo, the measure irf 
approaches the uniform measure ir^. For each value of the parameter t, the measure irf is 
invariant under the subgroup of rotations which leave invariant the vector ip, making it the 
preferred state of the measure. Our construction is based on the unitary Brownian motion, a 
stochastic process valued in the set of unitary matrices. At fixed time t, this process itself is an 
interpolation between the identity matrix (at t = 0) and the unique invariant Haar measure 
on the compact group of unitary matrices (when t — > oo). The construction is motivated by 
the similar procedure that was used in the definition of the Fubini-Study measure. 

The paper is organized as follows. In section 2 we review the definition and some basic 
properties of the unitary Brownian motion. Section 3 contains the definition of the new family 
of ensembles of pure states, that are analyzed in section 4 using the Laplace transform. Finally, 
we compute in section 5 averages of some quantities of interest in quantum information theory. 



Let us now introduce some notation. In quantum information theory, any norm one vector 
(or pure state) x gives rise to a probability vector. More precisely, if e$ is the canonical basis 
of H ^ C^, then we can decompose ip in the following form 

ip = y^^jej. 

To such a state we naturally associate the probability vector 
(1) P {iP) = (\iP 1 \\...,\iP N \ 2 ). 

Physically, if ej determines the level of energies of an atom and if x represents the wave function 
describing this atom, the quantity \tpi\ 2 represents the probability to be in the energy level ej, 
that is 

P[to be in the state ej] = iV'il 2 - 
Another physical motivation related to probability theory concerns the measurement of 
observables. It is known in quantum mechanics, that a physical quantity of a quantum 
system % ~ C N is represented by an observable, which is an auto-adjoint operator on %. 
Let A be an observable and A = Yli=i ^i^i be its spectral decomposition. If ip is a reference 
vector state of %, it follows from the axioms of quantum mechanics that a measurement of 
the observable A gives a random result Aj: 

P[to observe A;] = ||P^I| 2 - 

In particular if the projectors Pi are the one dimensional projectors on Ce^ we recover the 
previous probability. 
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Let us now recall some elements of probability theory that we are going to use. If X and 
Y are two independent real Gaussian random variables of mean and variance 1/2, then 
Z = X + iY is said to have a complex Gaussian distribution of mean and variance 1. We 
denote by A/"c(0, 1) the law of Z . A complex vector (Z±, . . . , Z n ) is said to have a multivariate 
complex Gaussian distribution A/j? (0, I n ) if the random variables Zi, . . . , Z n are independent 
and have distribution A/c(0, 1). 

We shall also extensively use the Haar (or uniform) measure Haarjy on the unitary group 
U(C); it is the unique probability measure which is invariant by left and right multiplication 
by unitary elements: 

(2) VV, W € Wjv(C), V/ : U N (C) -»■ C Borel, 

f f(U)dHa&r N (U) = [ f(VUW)dHaar N (U). 

JU N (C) Ju N (C) 



2. Unitary Brownian Motion 

This section is devoted to the presentation of the unitary Brownian motion: definition, 
properties, stochastic calculus, invariant measure. In particular, we present all the ingredients 
that we are going to use for generating random quantum states. 

The unitary Brownian motion refers to the natural definition of a Brownian motion on 
the unitary group of complex matrices. This is a special case of a Brownian motion on a 
differential manifold and more precisely on a compact Lie group (in differential geometry this 
is sometimes called the heat kernel measure). 

2.1. Definition. We now define the unitary Brownian motion (UBM). To start, let us intro- 
duce some notations. For N € N, we denote Un(C) the unitary group on Mn(C), that is, 
U N {C) = {U £ GL n (C) /UU* = 1} and let denote by M S ^(C) the set of Hermitian matrices 
on M N (C), that is, M%(C) = {H € M N (C)/H* = H}. The set A4^(C) is a real linear 
subspace of 7Wat(C), that we endow with the scalar product 

(A,B) = NTt[A*B] = NTr[AB]. 

This way, we can consider the Brownian motion on M S ^(C), the unique Gaussian process 
(H t ) which satisfies 

(3) Vs,t, Vi,Be M%(C), E[(A,H a ){B,H t )] = (s A t)(A, B). 

In an equivalent way, the process (Ht) has the same distribution of the random Hermitian 
matrix whose upper-diagonal coefficients are /\ N (Ht l + iCt 1 ) an d whose diagonal coefficients 
are — =Z?^, where (B^ 1 , C^ 1 , D%) are independent standard real Brownian motions. 

We have now all the ingredients to define the unitary Brownian motion. This is the process 
(Ut)t^o, solution of the stochastic differential equation 1 



(4) j uu t = i(dH t )U t - -U t dt 




1 

I. 



This equation is written in Ito form. 
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In particular, we have 

(5) du; = -iu* t dH t - ^u;at. 

As a warm-up calculation, let us check that the process (Ut) is unitary. In order to compute 
d(UfUt), we need to use Ito stochastic formulas for matrix valued stochastic processes. Such 
formulas have been derived in [1, Section 2.1]: if (X t ) and (Y t ) are two matrix valued stochastic 
processes defined by 

dX t = AtdHtBt + C t dt, dY t = DtdHtEt = F t dt, 
then we have the following Ito formula 

(6) d{X t Y t ) = X t dY t + dX t Y t + ^Tv[B t D t ]A t E t dt. 
This way, we obtain 

dUtUt = UtdUt + ^Ti[U t Ut}Idt 

(7) = (±Tr[U t U t *]I-U t U^dt. 

Since the unitarity condition is satisfied at t = (JJqUq = I) and the identity is a solution of 
the above ordinary differential condition, by uniqueness of solutions we have UtU£ = I for all 
time t. 

2.2. Laplace Beltrami operator, Markov generator. For the sake of completeness, let 
us describe the Markov generator of this process. In particular, this allows us to motivate the 
definition of the Brownian motion from a geometric point of view. To this end, we denote by 
utv(C) the Lie algebra of Un{C) and by (Xi, . . . , X N z) an orthonormal basis of u(N). For all 
smooth functions F : Un{C) — > R, we define 

(8) (£ Xi F)(U) = j^ =o (F(Ue tx >)) , Vi 
The operator ^A, where 



(9) a = ^£; 
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is then the Markov generator of the unitary Brownian motion (this justifies the name heat 
kernel measure which is sometimes used to defining this process). The operator A is actually 
the Laplace-Beltrami operator on the Riemmanian manifold Un{C) endowed with the Riem- 
manian metric induced by the scalar product on matrices {A, B) = N Tv[A* B]. Let us stress 
that this operator does not depend on any particular choice of an orthonormal basis. The 
Markov generator character of A is expressed in the following proposition. 

Proposition 2.1. [13, Proposition 2.1] Let F : Un{C) — s> C be a function of class C 2 . Then 
for all t ^ 0, we have 

(10) F(U t ) = F(I) + V / (£ Xt F)(U s )d(X k ,iH s ) + / -AF(U s )ds 

i=1 Jo Jo l 
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and the processes ((Xk, iH t )), k = 1, . . . , N 2 are independent standard real Brownian motions. 

The proof of this proposition relies on the classical Ito formula; this is a classical result of 
stochastic analysis on manifolds. A particular consequence of this proposition is that for all 
smooth function F, the process defined by 

F(U t ) - F(U ) - £ ^AF(U s )ds, 

for all t is a martingale (with respect to the natural filtration associated with the Brownian 
motions (B^ 1 ,C^ 1 , D^)) and Ut is the unique process "in distribution" satisfying such a prop- 
erty (this property characterizes the unitary Brownian motion and could be used as a starting 
definition) . 

2.3. Invariant measure. As announced in the introduction, the UBM will help us to define 
a new family of random states which interpolates between deterministic and uniformly dis- 
tributed random states. This relies on the large time behavior and the invariant measure of 
the UBM. 

Theorem 2.2. For all initial unitary conditions Uq the solution of the stochastic differential 
equation 

(11) dU t = i(dH t )U t -\u t dt, 

converges in distribution to the Haar measure Haar^ on Un(C). In other words, the Haar 
measure is the unique invariant measure of the Markov process, solution of (11). 

The theorem above justifies the property of interpolation between the identity operator 
(Uq = I) and the Haar measure for large time (t goes to infinity) for the unitary operator Ut- 
This property is essential for our definition of new ensembles of random pure states. This fact 
will be made precise in Section 3. 

Theorem 2.2 also shows that in general the distribution of the UBM (Ut) is not invariant 
by unitary multiplication (except under the invariant measure) but the distribution of (Ut) is 
nevertheless invariant by unitary conjugation and by inversion. 

Proposition 2.3. Let (Ut) be the UBM defined by the SDE (11) and let V be any unitary 
matrix in Un(C). The processes (VU t V*) and (Uf 1 ) have the same distribution than (Ut). 

Proof The property concerning the stochastic process (VUtV*) follows from the fact that it 
satisfies the same stochastic differential equation (11) with the same initial condition. Let 
(W t ) defined by W t = VU t V* for all t, we have 

dW t = V(i ■ dH t U t )V* - -VU t V*dt 

= i-d(VH t V*)Wt-^W t dt. 

Since (VH t V*) is a Brownian motion on the Hermitian matrices, we see that (Wt) and (Ut) 
satisfies the same SDE. Hence, as they start with the same initial condition, the two processes 
must have the same distribution. 

The statement for the inverse is an easy consequence of the Ito formula (11) and the fact 
that the inversion corresponds to the complex adjoint and then is a linear mapping (not 
affected by derivation) □ 
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2.4. Useful formulas. We continue by investigating some properties of the UBM which are 
going to be useful for studying the random pure states generated by the UBM. In particular, 
we will be interested in the properties of the coefficients of the matrix (Ut) which we denote 
by u( k ' , for 1 ^ j, k ^ n. The stochastic differential equations satisfied by these elements are 
of the following form 

N 1 

(12) dui k =ij2( dH h u t sk -^ui k dt. 

s=l 

and for their complex conjugates 

N -, 

(13) dU? = UfdH? - "-Utdt. 

s=l 

In the next section we will need the following expressions 
d\U^\ 2 = dUl 1 TlJ 1 + U^dUJ 1 + dUfduj 1 

(14) =iYXdH>')uM 1 -iJ2utU t sl dH^ + ( -\ut\ 2 + dt. 

s=l s=l ^ ' 

The previous formula relies on the stochastic bracket for the elements of (Ht), that is, 

(15) d((H?,H»)) = ^6 il 6jk, 

where S is the Kronecker delta symbol. From this, we obtain the brackets 2 of the matrix 
coordinates which is given by 

(16) ddpi^Au^)) = I (\u?\H jk - pi 1 | 2 |c/f | 2 )^. 

3. Random Pure States Generated by Unitary Brownian Motion 

Before developing our theory for random pure states generated by UBM, we review the 
definition and the basic properties of the uniform (or Fubini-Study) probability measure on 
the set of pure states (or unit vectors). The lack of any a priori information on the state tp 
of a quantum system described by a Hilbert space T~L ~ imposes the choice of a measure 
which should be invariant by changes of bases. In our setting of finite dimensional complex 
Hilbert spaces, changes of bases are implemented by unitary operators U G Un(C). As a 
consequence, we ask that the uniform probability measure should be unitarily invariant. A 
probability measure ir on the unit ball of T~i is called unitarily invariant if for all Borel subsets 
A and for all unitary operators U £ 

(17) it(UA) =tt(A). 

The above condition determines uniquely the measure ir: it is the normalized surface area 
of the unit ball of C^, which we shall denote by ir^, for consistency reasons which shall be 
clear later. Moreover, we introduce the image measure = sq_^7roo, where sq : C N — > An, 



2 Here we have adopted the notations (( , )} for the stochastic bracket not to be confused with the scalar 
product { , } . 
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sc l[(V'j)] = GV'il 2 )- I n other words, if the unit vector ip has distribution tt^, then the probability- 
vector {\tpi\ 2 )f = i has distribution a^. 

Other that the abstract definition of the invariant measure tt^, there are two more char- 
acterization of this probability that are important in what follows: 

(a) Let X € be a standard complex Gaussian vector. Then X/||X|| has distribution 

(b) Let U G Wjv(C) be a Haar-distributed random unitary matrix. The first column (or 
any column, or any line) of U has distribution -Koq. 

The first statement above is useful when one needs to sample from ir^. The second state- 
ment above will be the starting point for the definition of new probability measures on the 
unit sphere of C N . 

Start with a fixed vector ip € C N of norm one, and define the stochastic process (ipt), where 

(18) iPt = U t il>, Vt > 0, 

and (Ut) is a UBM starting at Uq = In- This gives rise to a stochastic process valued in 
the unit sphere, with ipo = ip. In the sequel, we study the properties of this process, whose 
distribution at time t we denote by trf. As before, the distribution of the probability vector 
\i/4\ 2 is denoted by o~f , making explicit the dependence in the initial condition ipQ = ip. 

Definition 3.1. The ensemble of pure states (unit vectors) of C N having distribution trf 
is called the unitary Brownian motion induced ensemble at time t. The distribution of the 
square moduli of the coordinates of a random vector ip in this ensemble will be denoted by 

Let us first discuss the connection between the distribution irf and the Haar measure tt^ . 
First, note that the distribution of ipt depends on the initial value ip, in contrast with Haar 
distributed random states, whose distribution is invariant. In particular, we do not have for 
nf invariance by all unitary transformations, that is, if V denotes a unitary operator, in 
general the processes (Vipt) and (ipt) have different distributions. More precisely we have the 
following result. 

Proposition 3.2. Consider a unitary operator V and a (possibly random) unit vector -ip. Let 
(ipt) be the process generated by a unitary Brownian motion independent of ip, with initial 
condition tpQ = ip. Then, the processes (Vipt) and (UfVip) have the same distribution. In 
particular, the processes (Vipt) and (ipt) have the same distribution if and only ifVip ~ ip. 

Proof. This follows from the fact that 

(19) Vip t = U t ip = (VU t V*)Vip = U t Vip. 

Since (Ut) and (Ut) = (VUtV*) have the same distribution and are independent of ip and 
Vip, the first part is then straightforward. The second part is a trivial consequence of first 
part. □ 

Corollary 3.3. Let ip be a random uniform unit vector, i.e. ip ~ ir^. Consider an inde- 
pendent UBM (Ut) which induces a process (ipt) with initial condition ip. Then, for all t, the 
vector ipt has distribution tToo, i.e. irf = -ir^. 
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Proposition 3.2 is more instructive when we look at deterministic initial conditions. In 
particular, the processes (Vipt) and (ipt) have the same distributions if and only if Vip = tp, 
restricting the class of unitary transformations which leave invariant the process. In other 
words, the distribution of (ipt) is invariant under all unitary transformations which fix the 
initial condition. In conclusion, the process (ipt) has a Un-i(C) invariance group, whereas a 
uniform unit vector tp ~ i:^ has a Wtv(C) invariance group. 

4. PDE for the Laplace transform of of 

We are now in the position to develop our model in more detail, using as a main tool a 
partial differential equation satisfied by the Laplace transform of the amplitude vector. Let 
tp € a fixed unit vector and consider the process generated by a UBM (Ut) with Uq = I: 

(20) <p t = U t ip. 

Consider also a fixed unitary matrix V such that Ve\ = tp, where e\ = (1, 0, . . . , 0) is the first 
element of the canonical basis of C^. The one can write 

(21) = U t V ei = U t e u 

where (Ut) is another UBM starting at Uq = V. Note that the choice of the matrix V satisfying 
Vei = ip is not important, because of Proposition 3.2. In this way, the initial condition of the 
problem has been transfered into the UBM (Ut) and we have 

iPt = U t e 1 = (U { 1 ,U?\...,U t N1 ), 

which corresponds to the first column of the unitary Brownian motion Ut- As is was discussed 
in the Introduction, such a unit norm vector gives rise to a probability vector (\tp 3 t \ 2 )j = 
(\U( l \ 2 )j having distribution of. This random variable is compactly supported, hence its 
Laplace transform determines its distribution. Let us define the Laplace transform by (notice 
the positive sign in the exponential) 



(22) ip(X u ...,X N ;t) = Ee^ 1 ^ = E 



N 



for all A = (Ai,...,Ajv) €E and all t ^ 0. Since the random variable is bounded, the 
function A i-> ip(\; t) is complex analytic for each t. 
The partial derivatives of the function ip read: 



(23) d j ip(X;t) = d Xj ip(X;t)=E 



V 3h2 e {x,\v t ^) 



(24) d jk ip{\;t) = d XjXk ip(X-t) =E [HI' 1 \ 2 \U t kl \ 2 e 

The following theorem is the main result of this paper, establishing a partial differential 
equation for the Laplace transform p. In principle, it allows to recover <p and then, by Laplace 
inversion, the probability vector (|C// 1 | 2 )^ 1 . 

Theorem 4.1. The Laplace transform ip of the random vector (If// 1 ! 2 ), j = 1, . . . , N satisfies 
the following partial differential equation 

(25) dtp = + ( - - A, V A ^) - - (A, H(ip)X) , 



RANDOM PURE STATES VIA UBM 



where A 2 represents the vector (X 2 )j = X 2 and \7\ and H represent the gradient and the 
Hessian operators, i.e. 

[ H {v)]jk = djd k ip. 

Proof. Using the multivariate Ito formula for the function (x\, . . . , xn) e^ X ' x ^ applied to 

the multidimensional process d^ 1 ] 2 ), we obtain: 

(26) 

N N 

dexp(A, l^ 1 ! 2 ) = ^A,exp(A, p^dp^ 2 + - £ A j A fc exp(A, [U^M |t/f| 2 , \U t kl \ 2 )). 

j=l j,k=l 

Taking the expectation and using (14), (16), we obtain 



dEexp(A, If/, 1 ! 2 } =jr XjK [exp(A, p 1 ] 2 } (-\U? \ 2 + -j-)] dt 

1 N r 9 

+ l - ^ Xj\ k E exp(A,|C/ t 1 | 2 )^(|[/f| 2 ^-|C/f| 2 |C/f | 2 ) 



(27) 



j,k=i 



dt. 



Using formulas (23), (24), we obtain immediately the announced partial differential equa- 
tion (25) satisfied by ip. □ 

From the above PDE one can obtain, in principle, all the information about the distribution 
o~f of the random vector (I?// 1 ! 2 ). In the remainder of this section, we shall focus on the 
marg inals If// 1 1 2 (j fixed); covariances p^^pf 1 ] 2 and other statistical quantities will be 
investigated in the next section. 

In the case of marginals, we compute in the following proposition the Laplace transform 
of the square modulus of one coordinate, in terms of the Kummer confluent hypergeometric 
function whose definition we recall : 



(28) 



iFi(a;6; z) = ^ 



k=0 



(a) kz k 

(b) k k\ ' 



where (x) n is the Pochhammer symbol, (x) n = x(x + 1) • • • (x + n — 1). 

Proposition 4.2. Let ipj(X;t) = tp(0, . . . , 0, A, 0, . . . , 0; t) = Eexp(A|C// 1 | 2 ) be the Laplace 
transform of the j-th coordinate of the first column of U t . Then 



(29) 

with the notation 
(30) 



A 



d m = J^Vj + ( — - A ) dyfij - —dxxfj, 



N 



A 2 



N 



d t <Pj 



dtpj dipj d 2 ifj 



dt 
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Given an initial condition I^q 1 ) 2 = c ^ [0, 1], there exists a sequence (a n ) n ^o of real numbers 
(depending on c) such that 

oo 

(31) tpj^t) = Y j a n e- Ant \ n iF 1 (n + l;N + 2n;\), 

n=0 

where 

, . n(n — 1) 

(32) An = n + V - ; . 

Proof. Equation (29) follows from equation (25) of Theorem 4.1 by letting A& = for all 
k j. In the rest of the proof, we shall drop the index j, since the initial condition will be 
encoded into ip(X;0) = exp(A|C// 1 1 2 ). 

Using separation of variables, we look for solutions of the form ip(X; t) = f(\)g{t). Neither 
of / or g can be zero, hence we obtain 

, ^ </ A 2 /" fx 2 \ f A 

Note that the left hand side of the above equation depends only on t and the right-hand since 
depends only on A. This is impossible unless both are equal to a constant C, in which case 
g(t) = e ct (we can move the constant factor to /) and / satisfies the ordinary differential 
equation 

(34) -» r+ (*- X )t + (*-c)f-*. 

Writing / as a power series /(A) = ^ n>0 a n \ n , we get: 

(35) Ca = 

(36) a 1 (l + C) = ^a 

(37) a J WLpl + k + c\ =^o fc _i Vfc>2. 
These equations can be summarized as (we put a_i = 0) 

(38) a k (C + A k ) = ^a fc _i VA: > 0, 

If C £ {— Afc}fc^o ; it follows that / = 0, which is impossible. Hence, C = — A n for some 
n ^ 0. We can compute all the coefficients of the series expansion of / from the recurrence 
relations above: 

(39) a m = ^ m < n 

(40) a n is free 

(a-\\ (n + l) m (n + m\ 1 

We obtain the final expression for the Laplace transform: 

°° °° /„ i ™\ \n+m 

n ^-Y j a n e-^ t Y j ( n + m \ X 

n=0 m=0 



(42) ^O-E-^AI^ I'J-E^E ' ; )W^7, 
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The second sum in the above formula admits a more compact expression using the Kummer 
confluent hypergeometric function iF±: 

and thus 

oo 

(44) <f(X;t) =Y,a n e~ Ant X n iF 1 (n+l;N + 2n;X). 

n=0 

□ 

Remark 4.3. Note that in the limit t — > oo, only the n = term survives and we obtain 

oo . m 

(45) Iim ¥ »(A;t) = X; 



m=0 

which is the result for the Haar measure. This is consistent with [11], Lemma 4.2.4, where 
the m-th moment in the Haar case was shown to be (^^v™^ 1 )- 

It is interesting to note that equation (29) does not depend on the actual value of j. 
However, it does depend on the initial condition c. Next, we compute explicitly (p for particular 
values of the initial condition c = | CT^ 1 1 2 € [0, 1]. The values of the coefficients a n appearing 
in the proposition can be computed in principle from the following initial conditions: 

ip(X;0) = e Xc 
<p(0;t) = 1. 

As a first observation, note that he latter relation fixes the value of the constant coefficient 
in the series, ao = 1. The first condition translates to (p = n + m): 

(46) ib a 4 P ) nvJo ^ =~v v ^°- 

The above infinite triangular system of linear equations can be solved in principle and explicit 
formulas for the coefficients a n can be found. 

Analytical formulas can be obtained (and easily proved by induction) in two particular 
cases. For c = 0, one can show that the unique solution to the equations above are given by 

(47 > -- (Jf + n 1 -!). - Vn> °- 

Similarly, for c = 1, one has 

< 48 > -- „ l (ff + ~,-i). - v " >a 

We were not able to obtain analytical expressions for all the coefficients in other particular 
cases. We gather next the first six coefficients in the important case c = l/N: 

fA^ 1 n -1 + iV 2(-2 + iV)(-l + AT) 

(49) oo = l, a 1= 0, °2-- 2N2{1 + N y <* = 3N s {2 + N){4 + N) . 
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(-1 + N) (30 - 29iV + 5iV 2 ) (-2 + N)(-l + N) (84 - 79iV + 77V 2 ) 

(50) a 4 - ~ 8iV 4(3 + iv)(5 + iV)(6 + AT) ' ° 5 ~ 15iV 5 (4 + AT) (6 + N)(7 + JV)(8 + N) ' 

5. Properties of the measure af 

This section contains a list of results which address important statistical properties of 
probability vectors distributed along the measure af. Moments, as well as covariances, are 
shown to satisfy ordinary differential equations that are solvable, see Propositions 5.1 and 5.3. 
We also compute quantities relevant to quantum information theory, such as average values 
of observables in Lemma 5.4 and bounds for average Renyi entropies. 

From Proposition 4.2, it is easy to obtain the expression of the moments of the j'-th coor- 
dinate \ip{\ 2 - To this end, we define a family of maps y p : [0, oo) — > [0,1], y p {t) = E|[// 1 | 2p . 
The maps, indexed by positive integers p ^ 1 depend implicitly on the size parameter N. 
The dependence on the index j is encoded in the initial condition y p (0) = \tpi\ 2p = \U^ 1 \ 2p . 
Interchanging the derivation operator d p and the expectation E, we obtain 

(si) %>(*) = ax£(0;t)- 

The applications y p satisfies a particular system of ordinary differential equations. 

Proposition 5.1. With the convention that yo = 1> the applications y p satisfy the following 
ordinary differential equations on R + : 

p 2 

(52) y' p = -A p y p + jjVp-i 

The solution of the system (52) can be expressed in terms of a n and A n in the following way 

(53) y P (t) = y( P V Ar a \ e- A "*. 

Proof. Using Proposition 4.2, we can rewrite the formula (42) for ip(X; t) in the form 

and then, taking the p-th derivative of the expression above, we obtain 



(55) y p (t)=p\Y^ 



P \ an c - Ant 

p-nj (N + 2n) p _ n 



n=0 



The coefficients a n depend on the initial condition (see the previous section). □ 

From the explicit computations in the previous section, we can specify to the coefficients 
a n for the initial condition ip = (1, . . . , 0). Indeed, the moments of \U} 1 1 2 , are 

P ^\ (N~l)n e - knt 



(56) E[\Um= P i± ( P _ ) 

n=0 VP UJ 



n) (N + n- l) n (N + 2n) p - 
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and the moments for l^ 1 ] 2 , j > 1, are given 



(57) 



E 



Pi 



l|2p 



^ Q \P - n) (N + n - l) n (N + 2n) p _ n 



-A„t 



For general initial conditions, one has to compute recursively all the moments y p . In the next 
proposition we give the general form of y\ and y2 for general initial conditions. 



Proposition 5.2. The moments y± and y2 are given by 
(58) yi{t) = 



V2(t) = 

(59) 

where yi (0) = \U^\ 2 



2/ 2 (0) 



1 



iV + 2 



iV + 2 



yi(o) 



AT/ N(N + 1) 



-{2+2/ N)t 



N 



e- L + 



N(N + 1) ' 



l^'l 2 an<iy 2 (0) = |l/f| 4 = |^'| 4 . 



We now move on to study the covariance of elements \U{ 1 \ 2 . We aim to compute the 
quantities 



E 



l|2|rrfcl|2 



, j?k = l,...,N. 

Note that for j = k, one can use the more general moment formula in Proposition 5.1. Let us 
define define a family of covariance maps fjk : [0, oo) — > [0, 1], with fjk(t) = E lU^l 2 ]^ 1 ] 2 
Using the Laplace transform, we have 



(60) 



We put also fj(t) = E lU^ 1 ] 2 ,Vi ^ 0. From Proposition 5.1, we have that 



/,-(*) = 0^(0;*) =E |C// 



)'1|2 



iV 



Proposition 5.3. For all j ^ k, the covariance applications satisfy the following first order 
ODE 



(61) 

which admits the solution 
fjk(t) = 

(62) 



fjk = Jf (fj + fk) ~ ( 2 + ^) hk. 



/i*(0)-]^2(/i(0) + /*(0)-^) 



N(N + 1) 



-t(2+2/JV) 



iV + 2 



/i(0) + /*(0) 



iV 



JV(JV + 1) 



£(0) = \^\ 2 , A(0) = |^T and f jk (0) = |^| 2 |VT- 
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Proof. To compute the two times derivative it is sufficient to consider the functions 

ip(\j, X k ;t) = ip(0, . . . , 0, Xj, 0, . . . , 0, Afc, 0, . . . , 0; t), 

with Xj (resp. Afc) in the j'-th (resp. fc-th) position. Indeed, we have fjkit) = d\ :j \ k ip(0,0;t). 
Using the Laplace transform of Theorem 4.1, we obtain 

(63) -i(Aja AjAjV + \ 2 t d, M f + 2A,A»8 AjArf ) . 

Applying d\ j \ k on both side and taking next Xj = Xk = 0, we get the expression 

(64) /j fc (t) = l(s AfcV (0,0;t) + a A . V (0,0;t)) - (2 + |) / jfc (t), 

which corresponds exactly to the expression (61). □ 

Naturally this expression depends on the initial condition. For example, in the case of 
tp = e\, we get 

hj{t) = (~¥+2 i 1 ~ n) ~ N(N + 1)) e ~ t(2+2/N) 



fjk(t) 



- ' » ~t{2+2/N) 



(N + 2)N N(N + 1) 



2 1 

(65) — : — e~* H ; -, for j 4 k and 7, k 4 1. 

v ' (N + 2)N N(N + l) J r J ^ 

Finally, we estimate different statistics of interest in Quantum Information Theory. More 
precisely, we compute the average value of an observable and we give a bound on the average 
Renyi entropy of a pure with distribution ivf . 

As stated in the introduction, the main motivation for this work was to define a new 
ensemble of random pure states. If a quantum system is in a state described by the vector ipt 
having distribution irf and an observable A £ _Mtv(C) is measured, quantum theory predicts 
that the average value observed is (ipt, Aif) t ) ■ 

We shall compute the average value E{ip t , Aip t ) = E{ip, U% AU t ip) = UtAU^ip). 

Lemma 5.4. The average value of the measure of a fixed observable A € Mn{C) on a 
quantum system described by the ensemble irf is 

(66) E<^,A^)=f(^^)-^P > ) e - t + ll ' ( ' l) 



N 



Proof. The result follows from the computation of dUfAUf with the Ito calculus. Using the 
Ito formula (6), we have 

(67) dU t AU; = U t {dAU;) + dU t (AU t *) + ^j^-Idt. 
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Computing d(ip, UtAU^ip) and taking the expectation, we get 
(68) 

which implies equation (66). 
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dE(^,U t AU^) = (-Ety,U t AU;i>) + dt, 



□ 



In particular, when t goes to infinity, we recover the usual result for the Haar measure 
7TQO. It is easy to compute average value of an observable A € Mn(C) for the Fubini-Study 
ensemble: 



(69) 
(70) 
(71) 



(A) = J (V-I^l^dvroo^) = J (ei|C/MC/|ei}dHaar([/) 

= <e x | J C/MC/dHaar(C/)|ei) 



TrA TrA 
= (e 1 \—I N \e 1 ) = —. 



We finish by deriving bounds for the Reyni entropies. Recall that these quantities are 
defined for a pure state ip by 



(72) 



N 

£i^i 2p 

3=1 



where p ^ 2 is an integer. For the random pure state generated by a unitary Brownian 
motion, we have 



(73) 



1 



S p (ip t ) = —\n 



P 



N 

i=i 



We shall estimate K[S p (^t)j with the help of the Jensen inequality. We have indeed 

N 



(74) 



1 



E[S P W*)] > ~, In 

1 — p 



1 

T~p 



E E i f/ * 1 i 2p 

J=l 

where Y p (t) = J2f=i ^\U^\ 2p are the sum-of-moments functions. Using the moment formulas 
obtained in the beginning of this section, one can compute in principle the function Y p (t). In 
particular, in the case when ip = ei, we have 

v 

(75) 



vm _ ,W P \ (N-l) n + (N -!)(-!) 



-A„t 



p-nJ(N + n- l) n {N + 2n) p _ n 
In the limit i — >■ oo, only the term n = survives, and we obtain 

lim F p (t) = (iV + p)f 



t— >oo 



iV + p 



which is formula (7.60) in [2]. We use this expression to bound the entropy of the coordinate 
vector of an uniform point on the unit sphere of C N : 



log 



1 (N+p\ 
N+p \ p ) 



p-1 
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Taking the formal limit p — > 1 in the above expression, we obtain the known bound [2, 
equation 7.70] 

N 1 

k=2 
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